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Abstract 

We consider directly the equations by which matter imposes anisotropies on freely propagat- 
ing background radiation, leading to a new way of using anisotropy measurements to limit the 
deviations of the Universe from a Friedmann-Robertson- Walker (FRW) geometry. This approach 
is complementary to the usual Sachs- Wolfe approach: the limits obtained are not as detailed, but 
they are more model-independent. We also give new results about combined matter-radiation 
perturbations in an almost-FRW universe, and a new exact solution of the hnearised equations. 
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1 Introduction 



The Cosmic Background Radiation (CBR) is the keystone of modern cosmological analysis, in par- 
ticular through use of the results of COBE [1,2] and other [3,4] measurements of anisotropies in the 
CBR to help us understand the nature of inhomogcncitics in the universe (sec [5] for a review of those 
observations). This paper presents a way of analysing this relationship that is an alternative to the 
usual analyses based largely on the Sachs- Wolfe effect (modified by astrophysical effects). Our analysis 
proceeds from slightly different assumptions than those usually made (though essentially compatible 
with them); it is to a considerable degree more model-independent than they are. 

In [6] we established a theoretical framework for investigating the direct implications of CBR 
anisotropies (i.e. those that follow without assuming particular inflationary or other evolutionary 
models for the universe). In that paper, we set up fully covariant and gauge invariant evolution 
and constraint equations governing the perturbations of the photon distribution and metric. These 
equations were then used to show that if 

Al: all fundamental observers measure the CBR to be almost isotropic in some domain, 

then it follows that 

A2: the spacetimc geometry is almost Friedmann-Robcrtson- Walker (TRW) - i.e. the 
shear, vorticity, spatial gradients and Weyl tensor are almost zero, and the metric may be 
put into perturbed-FRW form in that domain. 

The latter is the assumption which underlies the usual Sachs- Wolfe analyses - for it is the starting 
point used to set up the Sachs- Wolfe equations. 

In this paper, we use the formalism of [6], extending its results to find quantitative limits on 
the anisotropy and inhomogeneity of spacetime set by anisotropies in the CBR, without assuming a 
particular model for the origin of such perturbations. We end up with a series of estimates (Section 
4) relating the inhomogeneity and anisotropy of the universe directly to the background radiation 
anisotropies. In addition, we find some new exact results on perturbations in almost-FRW universes 
with both matter and radiation (Sections 2 and 5). For a class of almost flat-FRW universes, we re- 
duce the full set of linearised dynamical equations to two linear ODE's, and give a new exact solution 
at late times (Section 5). 

The difference from the more usual analyses is that they consider observations of the CBR made 
at one space time event P ('here and now'), relating them to assumed perturbations of the universe 
at the time of decoupling, these in turn taken to arise from some particular evolutionary history or 
other. Here we make no such evolutionary assumptions; however (cf. Al) we assume the nature of 
CBR anisotropies is known not only at our own space-time position P but also throughout an open 
domain containing that event (eventually chosen to represent the period between last scattering and 
the present, in a region containing our world- line and past light cone). 

At first glance this seems to make our analysis far more dependent on unobservable data than the 
standard approach. However this is an illusion, for that approach builds in equivalent assumptions 
at the outset, but in a rather hidden way, because it assumes A2, which cannot be proved on the 
basis of observational evidence alone [7,8]. It can only be deduced from such data on the basis of 
some kind of Copernican assumption such as Al, which is not directly testable [6-8]. Our approach 
helps make clciar precisely what these hidden assumptions are, and thereby enables one to test how 
weak they can be and still allow deduction of the desired results. In fact, as pointed out by George 
Smoot (private communication), the Copernican assumption is in principle partially testable. This 
follows since the Sunyaev-Zeldovich effect allows us to confirm to some degree that distant galaxies 
do see nearly isotropic CBR, for otherwise the scattered radiation would have a significantly distorted 
black-body spectrum. 
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In the covariant approach we work in the real non-FRW spacetimc, rather than starting from an 
assumed FRW model and perturbing away from it. Tlie former approach considers the full set of 
dynamical equations that govern the real variables. The latter approach risks missing certain effects 
[9,10] or masking underlying assumptions. We provide an example of the second kind by showing that 
the Boltzmann equation imposes constraints on the photon distribution if the monopole moment is 
assumed to be Planckian to first order. The covariant formalism for analysing fluid inhomogcneities 
[9,11,12] is a development of Hawking's approach [13], and gives a gauge invariant alternative to 
Bardeen's formalism [14]. The covariant approach to the photon distribution function in this paper 
and [6] is based on [15,16], and is an alternative to the application of a Bardeen-type formalism, as 
presented in [17,18] (in different contexts from that of this paper). 

Despite the success of the standard inflationary models with dark matter and critical density 
[19,20], current CBR observations are consistent with alternative models, and do not by themselves 
give independent tests of inflation [21,22]. Furthermore, future observations could produce problems 
for the standard models. The covariant approach provides a clear and direct relation between obser- 
vational and theoretical quantities, unobscured by particular gauge choices or by the complexities of 
harmonically determined variables. Furthermore, we do not impose any specific model to generate 
fluctuations in the CBR. Thus we investigate, as far as possible, what is determined directly by ob- 
servations of the CBR made by the family of fundamental observers. Where we are forced to make 
additional assumptions, they are made about observational quantities - and are thus in principle fal- 
sifiable by observation, provided we make some kind of Copernican assumption such as Al, stating 
that all fundamental observers see the same kind of thing (the nature of the required assumptions is 
clarified below). In this sense, we provide^ a framework for comparing and testing various models, in 
which there is as clear as possible a distinction between observed and assumed properties. 

Notation: 

The metric gab has signature (— , -h, -h, +)• Einstein's gravitational constant, the speed of light in 
vacuum, and Planck's constant are one. Round brackets on indices denote symmetrisation, square 
brackets anti-symmetrisation. Va is the covariant derivative defined by Qab- Given a four- velocity u", 
the associated projection tensor is hat = gab + UaUb, and the comoving time derivative and spatial 
gradient are 

Qa...b = U^VcQa-.-b , 
VcQa...6 = h/ha^ . . . hb^VaQe-.f 

for any tensor Qa...b (in [6] we used ^Va for Va). If the tensor is spatial, we define 

\Qa...b \={Qa...bQ'' -'Y'^ . 

Given a smallness parameter e, 0[N] denotes O(e^) and A ^ B means A — B = 0[2] (i.e. these 
variables are equivalent to 0[1])- When A ~ we shall regard A as vanishing (for it is zero to the 
accuracy of the first-order calculations that are the concern of this paper). 

2 Covariant and Gauge-Invariant Analysis 

The fundamental observers are assumed comoving with the cosmological matter, which is modelled 

by dust with mass-energy density p, and which is non-interacting with the radiation (as we are con- 
sidering the epoch after last scattering). The physically preferred four-velocity u" of this matter is a 
suitable average over peculiar velocities (which are small). The matter flow is characterised by and 
its rate of expansion O (= 3i? = 3S/S > 0, where H is the Hubble parameter and S the scale factor), 
shear aab, and vorticity ujab, all non-zero in general; however the flow lines are geodesic: = 
(consequent on the vanishing of the matter pressure). 

The frame defined by it" defines an invariant 3-1-1 splitting of tensors [23]. In particular, for a 
photon four-momentum 

p» = ^(w« + e»), e„u» = 0, e„e» = l, (1) 
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where E is the photon energy and Ca the direction of photon momentum, relative to fundamental 
observers. Then the CBR distribution function may be expanded as [15] 

f{x^, E, e") = F{x\ E) + Fa{x-, E)e'' + E)e''e'> + ... (2) 

where the covariant harmonics (multipole moments) Fai...ai (a;*^, -E) for L > 1 are symmetric trace- 
free tensors orthogonal to u", that provide a measure of the deviation of / from exact isotropy (as 
measured by u"). If the CBR is almost isotropic after last scattering for all fundamental observers, 
then [6]: 

F, F = O[0], VbFa,...a^=0[l\, L>1. (3) 

Energy integrals of the first three harmonics define the radiation energy density, energy flux and 
anisotropic stress [15]: 

POO 

11 = At: E^FdE = 3p = O[0] , (4) 
Jo 

Qa = Yj^ E^FadE = 0[1] , (5) 

TTab = ^j^ E^F,,dE = 0[l] (6) 
(in [6] we used jiR for ji). We will also need the integral of the third harmonic: 

Uc = ^ / E^FabcdE = 0[1] . (7) 
35 Jo 

Note that the total energy-momentum tensor is made up of matter and radiation contributions; 

Tab = (P + IJ')UaUb + ^l^hab + T^ab + "^U^aQb) ■ 

2.1 Covciriant linearised hcirmonics of the Boltzmann 
equation 

The Boltzmann equation in curved spacetime 

may be decomposed into covariant harmonic equations via (1), (2). The full (exact and non-linear) 
results arc given in [16, p501]. For the collision-free and zero acceleration case, the linearised zero, 
first and second harmonic equations are: 



BF 

EF-lE^e— + lEVaF''c^O, (8) 



BF 

EFa-\E^Q^+EVaF + lEybF\^Q, (9) 
EFab - \E^e^-^ - E^aJZ + E^iaFb)- (10) 



These are the fundamental (covariant) equations governing the dynamics of radiation at a micro- 
scopic level. 
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2.2 Covariant linearised evolution and constraint 
equations 

If (8-10) are multiplied by E'^ and integrated over all photon energies, then they produce the radiation 
conservation equations governing /i and qa , as well as the crucial evolution equation for TTaf, (given for 
the first time in [6], in full non-linear form). These equations and the remaining (linearised) conser- 
vation, Einstein, Ricci and Bianchi equations are as follows, obtained from [6] and the general exact 
equations of [23] (with Mo = but allowing for an imperfect energy- momentum tensor). 

a) Matter and radiation energy and momentum conservation: 

p+ep=0, Ua = 0, (11) 

ii+^^li + Vaq" -Q, (12) 

qa + |eg„ + \Va^l + ^b1^\ ^ . (13) 

b) Evolution of radiation anisotropic stress tensor: 

i^ab + |07ra6 + j^H^ab + 2V (aQb) - ^hab^ cQ" + ^cCab ^ . (14) 

c) Einstein, Ricci and Bianchi propagation equations: 

e + ie^ + ^u + ip-o, (15) 

O'ab + ^Qf^ab + Eab - ^TTab — , (16) 
Wab+feWab^O, (17) 

Kb + QEab + V''i/(„^£6)ed + (|p + |M)a„6 + (18) 

+ ^T^ab + 5©7ra6 + {aQb) " ^hab^ cq"" ^ , 

Hab + OHab - V'^£(/e6)ed + ^V^a'e^cd ^ . (19) 

d) Einstein, Ricci and Bianchi constraint equations: 

<?a-iV,e + V^((7b„+a;6a) = 0, (20) 
V„a;'^ = 0, (21) 

^fa6 + V''h„" + W(al£6M-0, (22) 

%E''a - iV„(M + p) + ^VfeTT^ + leqa ^ , (23) 

VbH\ -{P+ |M)Wa + i^abcV^^ ~ , (24) 

where Eab, Hab are the electric and magnetic parts of the Weyl tensor, and Sabc = Vabcdu''' , with rjabcd 
the spacetime permutation tensor. 

It follows [6] from these equations and (3-7) that: 

^, V- - O[0] , VaV- = 0[1] , C7ab, COab, Eab, Hab = 0[l] . (25) 

where ip = fj,, G, p. These qualitative results from [6] will be made more quantitative in section 4. 
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2.3 Integrability conditions and conserved quantities 

In the case of zero acceleration, the linearised form of the identity [11] governing the commutation of 
time and spatial derivatives is: 

~ {Van + ieVa^- (26) 

where ^' is any 0[1] spatial tensor or any scalar with 0[1] gradient. The commutation of spatial 
derivatives themselves is given by the projected Ricci identities [11], which imply the exact identity 

V[aV6]V = -IpLJab (27) 

where ip is any scalar, and for a nearly-FRW spacetime, the linearised conditions 

^[a'^b]1pc ^ ^K[atpb] (28) 

where tpa is any 0[1] spatial vector, and fc = 0, 1, —1 is the spatial curvature index of the limiting 
(background) spacetime; and 

V[aV6]V'cd ^ -^{hc[a'>Pb]d " V'c[a/i6]d) (29) 

where tpab is any 0[1] spatial tensor. These identities could be overlooked in an approach that starts 
from an FRW background solution and perturbs away from it. The integrability conditions implicit 
in (27 29) lead to the new results: 

(I) if the covariant vector perturbations are spatially homogeneous to first order, then the 
vorticity vanishes to first order (i.e. non-zero terms are at most second order), and either 
all vector perturbations vanish to first order, or the spacetime has a flat FRW background. 

(II) if the covariant tensor perturbations are spatially homogeneous to first order, then 
the spacetime is either FRW to first order, or it has a flat FRW background. 

The first result follows from (27), which shows that Uab — since Va^/" is a vector perturbation 
for ip = fi,p, and from (28), which implies kS~'^ipa — for ipa = Qa or any other vector perturbation. 
The second result follows from (28) and (29), which imply 

kS-^i>a ^ ~ kS-^^ab 

for any vector perturbation ipa (since VaV'b is a tensor) or tensor perturbation tpab- Thus either 
k/S'^ ~ (the background FRW model is flat to the accuracy we are working), or 

In the latter case, (27) implies uiab — 0, (14) implies aab — 0, and then (20) gives Va6 — 0, while (16) 
and (22) give Eab — ~ Hab. Thus the spacetime is FRW to first order (i.e. it differs from FRW at 
most by second order terms). □ 

We can derive further linearised integrability results that hold in general (i.e. without assuming 
homogeneity of vector or tensor perturbations), by covariant differentiation of the dynamical equations. 
For example, taking the gradient of (24), and using (21) and V^a^bQc] — (which follows from (28)), 
we get 

V„V6if"^~0. (30) 
Similarly, using the contraction of (29) for '^cd = i^cd, the gradient of (20) gives 

V„V6a«^ + Vaq" - i V^e ~ , (31) 

while (23) yields 

V„ V6 {E"-^ + ^tt"*") + i e Vag" - i V' (m + p) ~ . (32) 
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Finally, we note the existence of various quantities that are conserved to first order along the 
matter flow. For example, (17) immediately shows that 

while (26), (19) and (16) imply that 

Hab + ^^ab — ^' Hab = Hab + 'V'^Cr{a''Sb)cd- (33) 

Then (33) and (11) give 

It therefore follows that if either the vorticity or H*^/p are negligible (i.e. 0[2]) at last scattering, 
they remain so at all subsequent times. 



3 Temperature Anisotropy of the CBR 

It is important to realise that the covariant dipole moment Fa of the CBR distribution (sec (2)), 
although dependent upon the choice of cannot be set to zero by this choice, since it is frequency- 
dependent and its vanishing implies special conditions on the anisotropy of /. Since the u"-frame 
is physically defined by the matter, and is already assumed to have been corrected for local peculiar 
velocities. Fa represents a possible residual intrinsic frequency dependent dipole moment of the CBR 
distribution relative to the matter, with invariant significance. This distributional dipole however 
contains much more information than the dipole of CBR temperature anisotropy, which is in fact 
proportional to the energy flux qa given by (5) (see (41) below). 

Even for the temperature dipole however, one cannot separate the intrinsic dipole from that induced 
by pecuHar velocity of the observer [19,24]. It is standard to assume that the intrinsic temperature 
dipole is negligible after correction for peculiar velocity. This is equivalent to the non-trivial assump- 
tion that the average four-velocity of the matter coincides with the energy-frame [12] four-velocity 
of radiation. Although it can be justified for adiabatic perturbations within the standard model [24], 
we will not make this special assumption, so that we allow for an intrinsic dipole in the temperature 
(i.e. non-zero qa) after correction for local peculiar velocities. This approach accommodates future 
improvement of observational results for the peculiar velocity which arc independent of the CBR ob- 
servations, and which may reveal a non-negligible residual temperature dipole. From this point of 
view, the current limits on the intrinsic dipole should be related to the current uncertainties in the 
local peculiar velocity. 

Another aspect of the dipole moment Fa which appears not to have been previously recognised is its 
link to deviations from a thermal Planck spectrum in the monopolc moment F. This aspect is hidden if 
one starts from a background FRW solution and then perturbs - rather than considering the real non- 
FRW solution. The latter approach shows via the Boltzmann equation that non-trivial constraints 
are imposed on the dipole if F is Planckian to first order (as strongly indicated by observations [5]). 
For suppose that 



F(x",f;) ~2 



(34) 



where k is Boltzmann's constant. Then (34) and the Boltzmann monpole harmonic equation (8) imply 

^pa..(T^s\ 3E/kT 

-[t^ S) l-cosh(£/fcT)- ^'^^> 

Thus the dipole moment is subject to the restriction (35) if the monopole moment is Planckian to 0[1]. 
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One consequence of (35) is 

^ = -| + 0[l]. (36) 

In contrast to many other treatments (where TS = constant), T is not a fictitious background tem- 
perature, but is the gauge- invariant average temperature in the actual spacetime. Now observations 
of the CBR measure temperatures in different directions on the sky. The full-sky average temperature 
T{x°') at event x" is determined by the monopole harmonic of the photon distribution: 

/•oo 

liix") = a[T{x'')f = An E^F{x'',E)dE (37) 
Jo 

on using (4), where a is the Stefan-Boltzmann constant. A directional temperature is determined by 

all the harmonics (2) via the directional energy density per unit solid angle that is defined by the 
integrated (bolometric) brightness [25] (sec also [26-28]) 

/(a;«, e*-) = / E^f{x'', E, e^)dE = — [T(a;«) + (5T(a;«, e*)]^ (38) 
Jo 47r 

which defines the gauge-invariant fluctuation ST{x"', e^). 

The covariant multipole moments Tai...aL{^^) (-^ ^ 1) of temperature anisotropy are trace-free, 
symmetric tensors orthogonal to u", defined by 

T = ^ = r„e» + Tabe'^e" + Tabcc'^e'e' + ... (39) 

By (38), (37), (2) and (4) they are given in general, to a good approximation, as normalised integrals 
of the covariant distribution harmonics: 



'ai...aL — 



47r / E^Fa,...a,dE. (40) 



1 

These moments give a covariant and gauge - invariant description of the CBR temperature variation, 
with spectral information integrated out. By (5-7), (40) gives the dipole, quadrupole and octopole as: 

^ 3ga ^ ISTTofe ^ 35^a6c / . , x 

4/i Sfj, 8/i 

The 

Tai...aL a frame-independent alternative to the usual multipole coefficients Alm in an 

expansion in spherical harmonics Ylm- If we choose a standard triad in the rest space of u"' such that 
e° = (0,sin6'cos(/), sin6'sin0, cos^?), then the two approaches are linked (cf. [15]) by 

oo L oo 

The correlation function 

C{a) = {T{x,e)T{x,e')) , e"e^ =cosa, 

where (. . .) denotes a statistical average, is the key quantity in actual observations. In this paper we 
will not consider the detailed statistical analysis of the correlation function, which is given for example 
in [19-21,24,26,29 31], where an inflationary model for perturbations is assumed. Our concern here 
is with the underlying principles of how to relate observational limits to properties of the spacetime 
geometry in a model-independent way. 

Current CBR observations place Umits on Ta...b{to,y\c) where to is the proper time along our 
worldline C since last scattering and y\c are comoving coordinates of C. By Al, these limits may be 
extended to hold on each worldline y at a proper time to along the worldline after last scattering: 

I Tai...ai(to,y) |< eL{to) ■ 

As in [6], we assume that anisotropics are 0[1] back to last scattering, and extend the limit to hold 
for all times <t <to, thus obtaining the assumption 
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Bl: there exist 0[1] constants such that eL{t) < cl- Hence for any event after last 
scattering 

I Ta,...ai^ \< CL • (42) 

In principle (and possibly not too far off in practice), observations place limits on the comoving time 
derivatives of the multipoles. As before we assume 

B2: there exist 0[1] constants e^,e**,e*** such that 

I fa,...a, \<eie, \fa,...a,\ < e^O^ , |f„^...„J < 61**0^ , (43) 

where we have normalised the derivatives relative to the expansion of the universe (recall that > 
is O[0]), and we will not need the higher derivative limits. Since it is effectively impossible for us 
to move cosmological distances off C, there will not be direct observations of the spatial derivatives 
of the multipoles. However we will need limits on spatial and mixed derivatives up to third order, 
and so assume also, on general plausibility grounds (given the basic Copernican assumption), that 
additionally: 

B3: there exist 0[1] constants e^, e^, e^' such that 

^cTa,...a, I < , |VdVeTa,...aJ < e^O^ , | Ve VdVeTa,...a J < e'^!®^ , (44) 

and constants e'£ , e'£* ' such that 

|(VcT„,...„J-| < e2@' , |(Vcr„,...„J--| < ere' , |(VdVeT„,...„J-| < e^e^ . (45) 

One should note that we expect all the quantities e defined here to be very small: probably 
at most 10"'*. The exact isotropic case considered by Ehlers, Geren, and Sachs [32] corresponds to 
Cl = €l = c'l = 0, and is a special case of what follows. In fact, our special case is a small generalisation 
of the EGS result to include non-interacting dust matter in the source of the gravitational field. 



4 Model-independent Limits on Spacetime 
Anisotropy and Inhomogeneity 

The observationally based limits (42) on temperature anisotropy lead directly via (41) to the limits 

kal < ImCI , kabi < ^^£2 , \Lbc\ < ^/"CS (46) 

on the radiation anisotropy tensors. Limits on the derivatives of the radiation tensors arise from 
differentiating (41), using (42-45) and the evolution and constraint equations of Section 2.2. For 
example, the limits on derivatives of qa that we will need are 

\qa\ < l^edei + el) , \Vaqb\ < l^ee'i , (47) 

\qa\ < 4f f^e\20 + 40fl + 2f2M)ei + 2Ael + QeJ*] , |V„V69c| < i^^&^e'| , (48) 

|(V„g6)'| < lf^e\Se[* + 4e'i) , |(V„V6gc)'| < |^^eH^^T + 4e'i') , (49) 
with similar expressions for nab,^abc- The density parameters are 

^R = ^' = ^ . 

Now (46) and (47-49) (with their counterparts for iTab and ^abc) are used in the evoluiton and con- 
straint equations and their derivatives in order to find limits on the geometric and dynamic quantities 
that characterise the deviation of the universe from FRW form. We emphasise that the following limits 
are gauge-invariant, covariant and imposed directly by observational quantities, i.e. the temperature 
dipole, quadrupole, octopole and their derivatives. 
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From (13) and (14) we get immediately 

|V„m| _ ^ ^ ^^g^^ ^ ^2et + f 4) , (50) 



M < |e, + + 5e[ + f 4 . (51) 

The remaining limits require more complicated manipulation of the evolution and constraint equations. 
From (27), using V6(13) to get |VaV6/i|, we find 



^<9e;+3er + fe^ (52) 



From (23), using Vc(16), Vc(14) and [Vc(14)]' to get limits on \VcEab\, |Vc(Ta6| and |(VcCra6)'|, we 
find 

l^aPl . 27 rr,/ , f ^ 



Q < fHe'^ + j ff[12ei + 12et + 61e'^] + 

) F[165ei + 45ei* + llOca + GQe'^ + 9e2** + IScg] . (53) 
From Va(12), using (13)' and (26) to get |(VaAt)'|, we find 

< H[50ei + 51el + Qe*,* + 3e'/ + f + 18e^ + f e^*] + 4(20^ + r!M)i^ei . (54) 

From (16), using (14)" to get \&ab\ and (12), (15), (26), we find 

^ < H[50e[ + 15e'{ + f £2 + Me^ + 3^* + f + ^e'*] + ^(HOr + 15nM)ffe2 . (55) 

From (22), using Vc(27) (with ip = fj,) and VcV6(13) to get IVcWabI, and Vc(14) to get iVcUab], we 
find 



^ < H[A5e'[ + 9e'[* + 9e'^ + 3e'* + fe'^' + f e^"] . (56) 



These equations show explicitly the role of the dipole, quadrupole and octopole. For example, (50) 
shows that the gradient of radiation energy density or average temperature, which reflects inhomoge- 
ncous deviations from FRW spacctimc. is bounded by the limits on both the dipole and quadrupole 
of temperature anisotropy, while by (51), the shear, which reflects anisotropic deviations from FRW 
spacetime, is bounded by the limits on the dipole, quadrupole and octopole. 

Note that if we follow the usual assumption that the CBR temperature dipole is negligible (after 
correction for local peculiar velocities), then by (41), (42-45) it follows that all the ei's vanish: 

Ta-O ^ g„~0 ^ ei = e* = ... = e'/* = 0. 

In this case, the bounds in (50-56) are reduced, so that a negligible dipole reduces the limits on 
anisotropy and inhomogeneity. 



Using (12) and (37), (50) may be re-written 

|VaT| 



|T| 

from which it follows that 



< iei+3et + f 4, 



l^aTl \T\ 



T T 
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where 

T T 
dR = — , tR = ^ (57) 

IVaTl' |T| ^ ' 

are characteristic length and time scales defined by the CBR. In practice e^,e" ^■'^d especially e^, 
e'l, - ■ ■ are not known from observations. In order to produce more useful versions of the results, we 
need a reasonable estimate of these quantities. First, we make the reasonable assumption 

CI: the spatial gradients of the temperature multipoles are not greater than their time 
derivatives: 

/ ^ * // /* ^ ** /// /** //* ^ *** 

Next, wc can estimate e^, e^*, e^** by invoking the characteristic time tn defined by (57), leading 

to assumption 

C2: the bounds on the time derivatives of the temperature harmonics are estimated 
by e ~ eL/tR , Q^l* ~ el A| , O^e^* ^ ei/^ , so that, using (36) 

* 1 ** ^1, *** 1 , 

ti — , _ gtL , — 27^L ■ 

We can now use CI— 2 to re-cast (50-56) in terms of the observationally realistic eL- 



II T 



AL^<H{12e, + fe,), (58) 



e 



< fei + 3e2 + f 63 , (59) 



l^o-bl ,10, I 2 , |'f^r|^ 

-Q- < T^i + 15^2 , (6U) 



M ^ I77e2 + ^ f6n., + 1.34.0 + fi.ol + 1^ 

|Vae 



er < i^^^ + i ) [60^^ + 1^4^^ + 6^3] + ( ^ 1 i/[16ei + f e^] , (61) 



< Hm^ei + 862) + 4(2fifl + nM)Hei , (62) 

En. 



e "'^3 

'55, I 103, I 23, \ I 4 



e 



< if(f £1 + ^^^£2 + f ea) + ^(ll^fl + l5nM)H€2 , (63) 
^ < H{fe, + f|e2 + ^es) . (64) 



If the dipole is neglected, then we can set ei = in (58-64). Note also that for the nR,rtM and H 
on the right sides of (50-56) and (58-64) we may use the O[0] values, i.e. the values they take in the 
limiting (background) FRW spacetime, which has non-interacting dust and isotropic radiation - exact 
solutions are given in [33] . 

The bounds (58-64) are the main results of the quest for a direct link from feasible observational 
limits on the CBR to limits on the deviations of the universe from FRW since last scattering (within 
our past light cone). They are based on the reasonable assumptions Bl— 3, CI— 2. These results 
give for the first time a direct relation between CBR observational limits and limits to anisotropy 
and inhomogeneity, without assuming a specific evolutionary model, or the curvature index k of the 
background. 

From these limits we can obtain conservative estimates of present-time bounds on the anisotropy 
and inhomogeneity of the universe. Let 

e = max(ei,e2,e3) 
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denote the upper limit of currently observed anisotropy in the CBR temperature variation, and take 
{^r)o < 1- Then (58-64) imply 



IVaMi 



< 17Hoe 



Fab I 

e 



< 6e , 







l^abl 

e 



<4e . 







'^""'^ <[|(OM)o + 56]ifoe, KlOHoe, 



e 



e 



'^"^^'^ <[4(OM)o + 77]ifoe, ( ^ | < C{Q)Hoe 



e 



(65) 

(66) 
(67) 



where C(fi) = 5 + 200/(f7Aif)o- The latter is a relatively weak limit, consistent with inhomogeneities 
in the largo scale mattcir distribution. If we are willing to assume that (il,M)o — 1 today, we get a 
reasonably tight limit from (67). However, the observational evidence points towards a range of values 
between 0.1 and 0.3 as more plausible [34]. Including the lowest limits implied by nucleosynthesis, we 
can represent the range of possibilities by a table of values: 



{CIm)o 


0.02 


0.1 


0.3 


1 


C{fl) 


10005 


2005 


672 


205 



As we go back in time towards last scattering, whatever value it has today, CIm will rapidly ap- 
proach 1. 

We note that if the dipole is neglected, then we can set ei = in (58-64) to obtain the revised 
"dipole-free" estimates: 



< 5Hoe 



\Oab\ 

e 



< 4e 



e 



< e 



'^""'^ <[|(l^M)o + 38]i?oe, <4J?oe, 



^ / 



e 



<41,fi„). + 2|H„., (^] < 



140 



(f^M) 







Hoe , 



(65a) 

(66a) 
(67a) 



\ ' / 

which are considerably sharper. We recover the EGS result [32] on setting e = (exact isotropy 
implies an exact FRW solution). 



5 Almost Flat-FRW Solutions 

As already pointed out, we are unlikely to have direct observational information about the spatial 
gradients of the temperature multipoles. In section 4 we used the estimate CI that the spatial 
gradients are not greater than the time derivatives. Here we consider a more stringent, but apparently 
not unreasonable, assumption on the spatial gradients. In effect we extend to the spatial gradient 
bounds the assumption C2 already made on the time derivative bounds - i.e. we replace assumptions 
Cl-2 by 

D: the bounds on both time and spatial derivatives of the temperature multipoles are 
determined via the characteristic scales (57) of the CBR: 
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Qel^^ ^ el^ \eL , Oe'^ ^ ^ ^ ^ < + le,)eL , (68) 
tR an 

where wc have used (36) and (58). It immediately follows from (68) that = 0[2], and hence by 
(44) and (41) 

Va96~0, VaTTbc^O, V a^bcd ^ , ■■■ (69) 

Thus from an apparently reasonable assumption on the spatial gradients, wc arc led to the vanishing 
at first order of vector and tensor inhomogeneities in the radiation. The point is that the radiation 
time scale tR is a zero order quantity (it exists for exactly isotropic CBR), while the radiation length 
scale cIr is first order - it is only finite when there are inhomogeneities in the CBR. If inhomogeneities 
in the temperature fluctuation r are determined by the characteristic inhomogeneity scale, then, as 
shown by (68), they are negligible in comparison with anisotropies. This in turn leads to a very 
restrictive condition: 



if the radiation energy flux and anisotropic stress are homogeneous to first order, 
i.e. if (69) holds, which will follow if D is true, then 

either the spacetime is FRW to first order, or the spatial curvature vanishes to first 
order and so the background has a flat FRW geometry. 

This follows as a special case of the results (I), (II) derived in section 2.3. Note that it may 
also be derived (but with greater effort) by considering the integrability conditions of the constraint 
equations (20-24) in the case that (69) holds. 

This result may be viewed as providing an alternative motivation for the almost flat-FRW model 

of the universe. Furthermore, we are able to reduce the system of dynamical equations in this model 
to a pair of linear evolution equations, for the shear and energy flux. Using the 0[l] identity (26) with 
(69), the spatial gradients of (13) and (14) imply 

VaVb/U ~ ^ tOab^O , VaObc ^ , (70) 

where we have used (27). Then (70) reduces (22) and the spatial gradient of (16) to 

Hab^O, VaEbc ^ . (71) 
By (69-71), the system of dynamical equations (11-24) closes at 0[1] and reduces to a sub-system 

for aab, T^ab, Eab 

(^ab + fQcTab + Eab - ^T^ab — , (72) 

TTab + |67ra6 + ^/XCTafe ~ , (73) 

E^b + QEab + {\p + ln^ab - l^T^ab ^ , (74) 

and a sub-system for qa, VaV (V' = Mi ©): 

qa + |e5„ + iV„/x ~ , (75) 
2Vae - 3g„ ~ , (76) 

Va(p + /x)-©'?a^O. (77) 

These sub systems represent a decoupling of the anisotropy and inhomogeneity, since (72-74) con- 
tain no spatial gradients (because Hab, ^j^ab and the acceleration vanish to the order of the calculation, 
this solution belongs to the 'silent universe' class recently examined by Mataresse et al [35,36]). In 
the special case where we assume that the temperature dipole is negligible, (41) and (75-77) show 

qa^O ^ V„Ai ~ VaP V^e ~ 

which, together with (69-71), shows that the spacetime is homogeneous to first order. Thus in the 
dipole-free case when D is assumed, the spacetime is Bianchi I to the accuracy of the calculation. 
Conversely, 
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if there is inhomogeneity in the radiation and matter, and if the CBR characteristic 
length scale determines the CBR temperature inhom.ogeneity (i.e. if D holds), then the 
dipole cannot be negligible after correction for peculiar velocities. 

From now on wc will assume that the dipole is not negligible, i.e. Qa ^ 0[2]. In (72 77), the 
coefficients /x, p and O may be given their FRW zero-order forms, which are in fact the solutions of 
(11,12) (using (69)) and (15). In particular, these equations imply the 0[1] Priedmann equation [6] 



- 3(/3 + m) - 0. 



(78) 



Before we consider the decoupling and solving of (72-74) and (75-77), we give the limits that they 
imply, using (46,47) and (68,69): 



|Va/^l 



< 12i?ei , 



|V„p| 



< 16 



e 



< 20ijiJei 



|0"ab| 

e 



< 3e2. 



(79) 



Note that a bound on Eab does not follow directly. These results sharpen the bounds given by 
(58,59,61,62) (recalling that CI— 2 have been replaced by D). In fact, the bound on the matter 
inhomogeneity is drastically sharpened at late times, when fij^ ^ ^m- The source of this is the 
disappearance of the gradient of the electric Weyl tensor which produced the H/^m term of (61). 
This gradient, along with all tensor gradients (i.e. gradients of vectors and tensors), drops out by 
virtue of D. In other words: 

if only scalar inhom,ogeneities occur (i.e. all vectors and tensors are gradients of 
scalars), then the matter inhomogeneities are rapidly suppressed at late times. 

A potential problem arising from (79) is that the limits suggest the matter inhomogeneity is much 
less than the radiation inhomogeneity at late times. However, the upper limits do not force this to 
occur, they simply allow the possibility. Below we will give an example of a late-time solution where 
the matter inhomogeneity is greater than the radiation inhomogeneity. 



By taking u"Va derivatives and using (15) and (78), the sub-systems (69-71) and (72-74) may be 
decoupled: 



qa + 3e<Za + (f - 2p 



_ 2 
3 



^^^p)<yab 

p)qa - . 



(80) 



(81) 



In principle we can obtain the 0[1] solution for the almost flat-FRW model after last scattering as 
follows. The solutions of (11), (12) and (78) imply 



rS- 



mS- 



e~5-2[3(r + rn5)]^/2^ r ~ = m 



which allow us to reduce (80,81) to linear ODE's in S. We write 

(7„6^4'6^S(,)(5), 7=1,2,3 

g„~Sf)g(A)(5), A = 1,2 
Then E(/) are linearly independent solutions of 



with ~ ~ 



(A) 



2(8r + 9mS) 



_ S{r + mS) 



r(94r + 89mS') 
_ 5'2(r + mS')2 



(82) 

(83) 
(84) 

(85) 



(480r2 + 1006rm5 + 525m^S^) 
20S'3(r + mS)2 



E~0 
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and (5(A) are linearly independent solutions of 

"3(4r + 5mS') 



Q" 



2S{r + mS) 



Q' 



2(5r + 7mS) 
S'^{r + mS) 



0^0. 



(86) 



Given the solution S(S') of (85), we find from (72,73) that 



[3 (r + mS*)] 1/2 



Cab 



^06 



[3(r + m^)] 
35 



1/2 



dS, 



[3(r + m5)]i/2 



(87) 
(88) 



where Cab — 0. 



Similarly, given the solution Q{S) of (86), we find from (75-77) that 
V„/x ~ -[Sir + m5)]i/2i3('^)[Q'(^) + iq^^^j , 



(A) 



(89) 

(90) 
(91) 



VaP ^ -[3(r + m5)]V2B(A)[Q'(^^ + -Q^^^] . 

Note that because two time derivatives were needed in the decoupling that led to (80), there will be 
a consistency condition imposed on the integration constants ^^p, Cab via (74). 

Thus (78), (83,87,88) and (84,89-91) represent an exact solution of the linearised equations gov- 
erning the metric, GBR and matter after last scattering in a universe subject to the observational 
assum,ption D, and with anisotropic radiaiion. The explicit analytic form of the solutions depends on 
finding explicitly the solutions to the ODE's (85,86). We can provide explicit solutions for late times 
(i.e. long after last scattering) when 

/X < p ~ « 1 (92) 

using (78). By (92), (80) simplifies to 

^a5 + l^'^ab + §Q^&ab + ^^^(^ab ^ . (93) 
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Writing aab — Uab&"', where Uab — 0, we find from (93) that n — 3, 3 
at late times. Thus an acceptable solution to (93) is 



, 2. However n = i violates (51) 



(Tab ^ UabS^/^ + VabO^ , Uab ^ ^ Kb ■ 



(94) 



Using (92) and (94) in (72-74) we find that the consistency condition on the constants of integration 
gives Vab — and the late-time solution is 



fab — Uab&^^^ , T^ab — \UabQ^^^ — h^CF ab , Eab — jUabQ^^^ — il^ab ■ 



(95) 



Now we use (92) to simplify (81) to 



qa+iQqa+^e\a-0- 



(96) 



With Qa — Ja^^, Ja — 0, we get n = |, |, neither of which violates the limit (46). Thus a solution of 



(96) is 



(97) 
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and this leads via (75-77) to 

v„M~-ij„ei°/3, VaQ^^Qa, v„p:^|j„ei°/3 + i^„e"/3. (gg) 

Thus we have presented an exact and exphcit late-time solution (94, 95,97,98), which in particular 
provides a quantitative measure for the rate of approach towards the limiting background solution. 
This solution also allows us to tighten the limit in (79) on shear at late times, and to give a late-time 
limit on the Weyl tensor: by (46) and (95) we have 

\aab\ . /160fi\ \Ea 



Furthermore, the solution contains spacetimes in which the late-time radiation inhomogeneity is 
negligible in comparison with the matter inhomogeneity. This arises when Ja and Ka in (98) are 
chosen to ensure that 

He, (100) 

where the final inequality follows from (79), consistently with (97,98). The particular choice of Ja — 0, 
i.e. radiation inhomogeneity vanishing to first order at late times, will achieve (100), with 

V„,.^0, ^^<^(^)He,, (101) 



M 



by (97,98) and (46). 



6 Conclusion 

We have seen how a series of assumptions of increasing sharpness (incorporating the inevitable Copcr- 
nican supposition) leads to increasingly powerful deductions from the CBR anisotropy. As emphasized 
before, these limits are independent of detailed assumptions about the dynamical history of matter in 
the universe, and provide an alternative mode of analysis to the usual approaches. This analysis has 
the advantage of being both covariant and gauge-invariant [11,12]. It gives somewhat less information 
than the usual approaches based on the Sachs- Wolfe effect and its generalisations, precisely because it 
is more model independent; however this also means that its conclusions are more robust than those 
more standard analyses, as they do not depend so much on the assumptions of particular evolutionary 
models. In particular they are not dependent on whether or not inflation took place, and whether or 
not the density parameter Cl is near the critical value. 

The analysis proceeds through a series of increasingly restrictive Copernican assumptions about 

the nature of CBR anisotropics in an open neighbourhood about our world line, which is envisaged as 
including the observable region of the universe. Such assumptions are inevitable if we wish to justify 
the assumption of an almost-FRW model [7,8]. The qualitative assumption Al (Section 1) is sharp- 
ened to the quantitative assumptions Bl-3 (Section 2) leading to the limits (50-56), which can be 
sharpened a little if we assume that the matter and radiation frames coincide (i.e. if there is no CBR 
temperature dipole). Somewhat sharper assumptions CI— 2 (Section 4) give better limits, leading to 
our estimates (58-64) and (65-67), which arc further sharpened if the dipole can be neglected to give 
(65a,66a,67a). These equations make quantitative the results of [6] (the universe is almost FRW, see 
A2) and include as a special case the Ehlers-Geren-Sachs exact theorem [32] . They confirm previous 
estimates based on the CBR anisotropics that the shear and vorticity are at most about 10~^ of the 
expansion (on choosing e ~ 10~*, to concur with recent CBR anisotropy measurements). 

We regard these assumptions and limits as highly plausible, and believe they are useful not only in 
terms of the limits obtained but also in making quite explicit the kind of Copernican assumptions one 
has to make in order to extract information from the CBR anisotropics (such assumptions necessarily 
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underlie the standard Sachs- Wolfe type analyses, because these analyses assume A2 as their starting 
point, but they do so in a somewhat hidden way). More debatable are the stronger assumptions D 
of section 5, which seem on the face of it quite plausible but then lead to very restrictive conclusions: 
either the universe is FRW to first order (that is, its difference from a FRW geometry is at most sec- 
ond order) or it has a flat background FRW geometry (to the accuracy of our first-order calculation). 
When this is true we can get explicit solutions of the equations representing the combined matter and 
radiation system - but they only allow spatial inhomogeneity when the dipole term cannot be neglected. 

Thus some may wish to adopt these stronger assumptions, while others may feel the conclusions 
are too strong and therefore the assumptions should be questioned. We are open minded in this 
matter; the main point is that the analysis presented here makes quite clear the range of possible 
assumptions, and their consequences when we take the Einstein-Liouville equations (and consequent 
propagation and conservation equations) into account. 

Finally we believe this paper shows well the utility of the covariant harmonic approach to both 
perturbations and to kinetic theory. In particular in examining kinetic effects, it makes quite clear 
how only the first three harmonic terms in the distribution function explicitly enter the field equations; 
the anisotropics represented by higher harmonics can affect the geometry only by cascading [15]: that 
is, by inducing anisotropics in the lower order harmonics through divergence or gradient terms, as in 
equation (10), or through collisions [16,37]. 
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